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Abstract. In this paper we show that a quasi-exactly solvable (normalizable or 
periodic) one-dimensional Hamiltonian satisfying very mild conditions defines a fam- 
ily of weakly orthogonal polynomials which obey a three-term recursion relation. In 
particular, we prove that (normalizable) exactly- solvable one-dimensional systems are 
characterized by the fact that their associated polynomials satisfy a two-term recursion 
relation. We study the properties of the family of weakly orthogonal polynomials de- 
fined by an arbitrary one-dimensional quasi-exactly solvable Hamiltonian, showing in 
particular that its associated Stieltjes measure is supported on a finite set. From this 
we deduce that the corresponding moment problem is determined, and that the k-th 
moment grows like the k-th power of a constant as k tends to infinity. We also show 
that the moments satisfy a constant coefficient linear difference equation, and that this 
property actually characterizes weakly orthogonal polynomial systems. 
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1. Introduction 



In a recent paper, [1], C. M. Bender and G. V. Dunne introduced a remarkable family 
of orthogonal polynomials associated to the one-dimensional Hamiltonian 

ff = _ a , + (4 £ ^4 £ -3)_ (4s + 4J _ 2) ^ + x e i (1) 

where J is a positive integer and s is a real parameter. If ift E (x) denotes an eigenfunc- 
tion of H with energy E, the polynomials P^{E) in question are proportional to the 
coefficients in the expansion of e x / 4: x^~ 2s ip E {x) in powers of x 2 , namely 

These polynomials are easily shown to satisfy a three-term recursion relation, from 
which it follows, [2], that they are orthogonal with respect to a certain Stieltjes measure 
du(E) (u being a function of bounded variation): 

J P k (E)P l (E)ME) = 0, k^l. (2) 

The form of the coefficients of the recursion relation satisfied by the polynomial 
system {Pk(E)}^ =0 implies that this system has several remarkable properties. First 
of all, the norm of the polynomials P k with k > J vanishes. Thus, the polynomials 
P k form what is called a weakly orthogonal polynomial system, [2] . To be precise, we 
shall use from now on the term orthogonal polynomial system for a family of orthogonal 
polynomials {-Pfc(-^)}fcLo wrt h degP fc = k for all fc, and such that the norm of P k does 
not vanish for any k. 

Secondly, each P k with k > J factors into the product of Pj and another polyno- 
mial, i.e., 

Pj+m = PjQm, ™>0, 

where Q m has degree m. Finally, the J simple real zeros of Pj are eigenvalues of H 
whose corresponding eigenf unctions, being the product of the factor 

\i\x) = e 4 x 2 

times a polynomial in x 2 , are square-integrable. The existence of these exactly com- 
putable eigenf unctions and eigenvalues of H had been deduced before, [3], [4], from 
the fact that H is quasi- exactly solvable, meaning that it is an element of the en- 
veloping algebra of a certain realization of sl(2, M) in terms of first-order differential 
operators acting on a finite-dimensional subspace of the space of C°° functions (see 
the next section for more details). The above results strongly suggest, [5], that there 
is a connection between quasi-exactly solvable Hamiltonians and certain families of 
weakly orthogonal polynomials. In this paper, we show in detail that this is indeed 
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the case for all one-dimensional quasi-exactly solvable Hamiltonians, both normalizable 
and periodic, satisfying very general conditions. The paper is organized as follows. 

Using the results on quasi-exact solvability reviewed in Section 2, we explain in 
Section 3 how to construct the weakly orthogonal polynomial system associated to each 
of the normal forms of a one-dimensional quasi-exactly solvable Hamiltonian listed in 
[4], [6]. Like the polynomial system introduced in [1], this system always satisfies a 
three-term recursion relation, whose coefficients we explicitly compute. This allows us 
to prove that one-dimensional (normalizable) exactly solvable Hamiltonians are char- 
acterized by the fact that their associated polynomials satisfy a two-term recursion 
relation. In Section 4 we show that the polynomials associated to an arbitrary one- 
dimensional quasi-exactly solvable Hamiltonian enjoy properties completely akin to 
those listed above for the Hamiltonian (1). We also study in this section the proper- 
ties of the moment functional defined by the family of weakly orthogonal polynomials 
of a quasi-exactly solvable Hamiltonian, giving a rigorous proof of the fact that its 
associated Stieltjes measure is supported on a finite set, [5], so that the integral (2) 
reduces to a finite sum. From this we deduce that the associated (Hamburger or Stielt- 
jes) moment problem is determined, and that the k-th moment behaves like the k-th 
power of a constant for large k, illustrating this statement with an explicit example 
for the Hamiltonian (1). We also show that the moments satisfy a constant coefficient 
linear difference equation, a property which in fact characterizes weakly orthogonal 
polynomial systems. The paper ends (Section 5) with a brief review of these results, 
stressing the role played by weak orthogonality — as opposed to true orthogonality — in 
their derivation. 

2. Quasi-exactly Solvable Potentials 

For the reader's convenience, we present in this Section a summary of the major results 
in the theory of quasi-exactly solvable systems that we shall need in the sequel. A one- 
dimensional Schrodinger operator (or Hamiltonian) H = —d\ + V(x) is quasi-exactly 
solvable if there exists a finite-dimensional Lie algebra of first-order differential operators 

= Span{£ a (z)d x + rj a (x) | 1 < a < r} = Span{T a (x) | 1 < a < r} 

such that: 

i) g leaves invariant a finite-dimensional module of smooth functions 91 C C°°(IR), i.e., 
X-f e 91 for all / e 91 and all X e 0. In other words, admits a finite-dimensional 
representation in terms of smooth functions. 

ii) H is in the universal enveloping algebra of 0, i.e., H can be expressed as a poly- 
nomial in the generators T a , < a < r, of 0. 

A Lie algebra of first-order differential operators satisfying i) is called quasi-exactly 
solvable. A Hamiltonian H satisfying condition ii) above for an arbitrary (not neces- 
sarily quasi-exactly solvable) Lie algebra is said to be Lie- algebraic. 

If H is quasi-exactly solvable, it follows that the restriction of H to 91 is a finite- 
dimensional linear operator 91 — > 91, and therefore the eigenf unctions of H lying in 91 
and its corresponding eigenvalues can be exactly computed by purely algebraic methods 
(diagonalizing a square matrix of order dim 91). We shall refer to these eigenf unctions 
of H as lying in 91 as its algebraic eigenf unctions (although, of course, they need not be 
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algebraic functions in the technical sense of the word). The functions in 9T need not a 
priori satisfy any boundary conditions (like square-integrability, periodicity, vanishing 
at the endpoints, etc.) coming from the physics of the problem, whose mathematical 
purpose is to guarantee that H is a self-adjoint operator. If they do, then the restric- 
tion of H to 91 is self-adjoint, and therefore H has exactly dim 91 linearly independent 
algebraic eigenfunctions, whose corresponding dim 91 real eigenvalues (counting multi- 
plicities) are exactly (i.e., algebraically) computable. We shall say in this case that the 
quasi-exactly solvable potential H (or the potential V) is fully algebraic. See [4] and [6] 
for an in-depth discussion of fully algebraic potentials under the boundary condition of 
square-integrability on R. 

It can be shown (cf. [7]) that a quasi-exactly solvable Schrodinger operator H can 
be expressed as a polynomial of degree at most two in the generators T a , < a < r, 
of g. Moreover, a well known theorem, [3], [8], [7], asserts that every quasi-exactly 
solvable Lie algebra of first-order differential operators g is related by a (local) change 
of variable 

z = C(x) (3) 

and a gauge transformation with gauge factor /j,(z) > to (a subalgebra of) one of the 
Lie algebras g n = f) n © R, where f) n = Span{J™, J n , J™} « sl(2, R), 

J- = d z , J£ = zd z --, J^ = z 2 d z -nz (4) 

and n is a nonnegative integer. In other words, every element X(x) Eg is of the form 

1 



X(x) = (jl(z) ■ J(z) 



C(*) 



for some fixed n. This implies that the gauge Hamiltonian 

1 

H{z) 



(5) 



x=t-Hz) 



is also a polynomial of degree at most two in the generators J", i.e., (dropping the 
explicit n dependence in the generators J™) 

-^gauge = °ab J a J b + Yl °a J a + C ^ ( 6 ) 

a,b a 

for some real constants c^, c a , and c ab = c ba (the minus sign is for later convenience). 
The spectral problems of H and H gauge are related in an obvious way: indeed, from (5) 
it follows that if x( z ) ^ s an eigenfunction of H gauge with eigenvalue E then 

*l>(x) = rtz)x(z)\ z=c(x) (7) 

will be an eigenfunction of H with the same eigenvalue (not taking into account the 
boundary conditions). Since the Lie algebra g n admits as invariant module the space 
V n of real polynomials of degree at most n in z, if H is fully algebraic then H gauge has 
n + l linearly independent algebraic eigenfunctions lying in V n . Hence H has n + 1 
linearly independent algebraic eigenfunctions of the form (7), with \ £ V n a polynomial 
of degree at most n. 

From (4) and (6) it follows, [4], that the gauge Hamiltonian is of the form 
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"iW= P(z)% + {«(*) - (*)} 9, 

+ {fl-iQ'W + ^P"w}, (8) 

where P, Q and P are polynomials of degrees 4, 2 and 0, respectively, given by 

P(z) = c ++ z 4 + 2c +0 z 3 + c 00 z 2 + 2c _z + c , 

Q(z) = c + z 2 + c z + c_, 
n(n + 2) 

Note that, due to the Casimir relation 

1 77 

J 2 --(J+J_+J_J + ) = I (n + 2), 

we have set, without loss of generality, c_| = 0. There are also explicit formulas for 

the change of variables (3) and gauge factor fx(z) needed to put the differential operator 
(8) in Schrodinger form, cf. [4]. Indeed, assuming that P(z) > on an interval I then 
for z E I we have 

* =c ~ lw= /*7fc' ^'^-"^{fMM (l2) 

and 

, (13) 

z=C(x) 

where the primes denote derivatives with respect to z. 

The canonical form (8) of the quasi-exactly solvable Hamiltonian H is not unique, 
since there is a residual symmetry group preserving the Lie algebra f) n , given by the 
adjoint action on f) n of the Lie group of transformations generated by g n = t) n ©!L More 
precisely, the elements of Q n are the infinitesimal generators of the standard GL(2, R) 
action on the space V n , given by 

p(z)eV n ^p(w) = ( 1 w + 6)y(^^y ^ 5)gGL(2,R). (14) 

We shall denote, as is customary, by p n this (irreducible) multiplier representation of 
GL(2,R) on V n . Note that the action (14) is just the composition of the projective 
transformation 

aw + (3 
z = r 

'fW + o 
5 



(9) 
(10) 

(11) 



V(x) = -R+ 



-n(n + 2) (PP" - |P /2 ) - 3(n + 1) (QP' - 2PQ') + 3Q S 



12P 



and the gauge transformation with gauge factor p(w) = ('jw + 5) n . The adjoint action 
of GL(2,R) on () n induced by (14) is given by 

J(z) ^ J(w) = { 1W + 5) n -J ( aw + l3 \ . ( 7W + 6 )~ n . (15) 

\ 'yw + o J 

A straightforward calculation, [4], shows that the generators of f) n transform under 
the representation p 2 ,-i — where p ni = p n <g> det 1 , det : A \— > det A being the standard 
determinantal representation — independently of n. As a consequence of all this, the 
transformed differential operator 

£ gauge = + ST ■ H gauge (^|) • (TV + S)~ n (16) 

is still of the form (8), with P, Q and R replaced by appropriate polynomials P, Q and 
R of respective degrees 4, 2 and 0. It can be shown, cf. [4], that R = R and 

p (w) = (jie±^1p ( E2t±e) , q (w) = (jmpl Q ( ?£±f\ , (17) 

with 

a (3 



det 



7 5 



Hence the polynomials P, Q and R determining the differential operator H gauge trans- 
form under the representations p 4 _ 2 , Pi -\ an d Pq of GL(2,R). Furthermore, the 
algebraic eigenfunctions of H gauge clearly transform under the representation p n ; in- 
deed, if x( z ) is an eigenfunction of H gauge with eigenvalue E then it follows from (16) 
that 

X(w) = (-yw + S) n ■ X (^4) (18) 

is an eigenfunction of H gauge with the same eigenvalue. 

In [6] and [4], the form-invariance of the differential operator H gauge under the 
GL(2, R) action (16) described above was exploited to place H gauge in canonical form. 
Indeed, it can be shown that there are ten inequivalent real normal forms for a (nonzero) 
fourth-degree polynomial P [9] transforming under the representation p 4 _ 2 of GL(2, R), 
each of which leads to a canonical form for H„ a „ ao . Of these ten canonical forms, five 
correspond to normalizable Hamiltonians, whose algebraic eigenfunctions are square- 
integrable (provided the coefficients c ab and c a satisfy certain inequalities), and the 
remaining are associated to Hamiltonians with periodic potentials. The five normal 
forms associated to normalizable Hamiltonians, which are characterized by the fact 
that P has at least one multiple root on the real projective line MP, are given by 

1. Hz 2 + i), 

2. z/(z 2 -l), 

3. vz 2 , (19) 

4. z, 

5. 1, 
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where v > is a real parameter. For example, the quasi-exactly solvable potential 
discussed in [1] corresponds to the fourth normalizable canonical form P(z) = z. The 
remaining normal forms, corresponding to periodic potentials, are 

6. z 2 )(l- k 2 z 2 ), 

7. v(l-z 2 )(l- K 2 (l-z 2 )), 

8. v(l + z 2 )(l + (l-K 2 )z 2 ), (20) 

9. v(\ + z 2 f, 
10. v(\-z 2 ), 

where v > 0, < k < 1. 

3. The Recursion Relation 

Let H = —d 2 + V(x) be a quasi-exactly solvable Hamiltonian. From the previous 
section, we know that there is a change of variable (3) and gauge factor fi(z) > such 

that H(x) = n{z) ■ H gauge (z) ■ with # gauge given by (6) (and c + _ = 0). 

Furthermore, if H is fully algebraic then it has n + 1 algebraic eigenfunctions of the 
form (7), with x( z ) e an eigenfunction of H g£iuge . Let Xe( z ) an eigenfunction of 
H gauge with eigenvalue E (not necessarily a polynomial in z). Writing 

oo 

XE (z) = Y / P k (E) Xk (z), (21) 



where 



fc=0 



z k 



and taking into account that 

J_ -Xfc = Xfc-i, ^o-Xfc=(^-|)x^ ^+ -Xfc = l)x fc+ i, (22) 

cf . (4) , we easily find that the coefficients P k (E) satisfy the following five-term recursion 
relation: 



-c__P k+2 = [(2k -n + l)c _ + c_] P k+1 

+ 



£? + C „ + c (fe--) +c 00 (A;--) 



+ - 1 - n) [(2k -n- l)c +0 + c + ] P fc _! 
+ fc(A; - l)(/e - 1 - - 2 - n) c ++ P k _ 2 , 



k > 0. 



(23) 



If c 0, the general solution of the recursion relation (23) depends on the two 

arbitrary functions P Q (E) and P ± (E). This simply reflects the fact that when c__ 7^ 
the leading coefficient P(z) of H gauge does not vanish at z = (cf. (9)); thus, the 
differential equation (H gauge — E) X e = nas a regular point at the origin, and therefore 
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it admits two linearly independent solutions (21) analytic at 0. If P (E) and P 1 (E) 
are chosen to be polynomials in E, then (23) implies that all the coefficients P k (E) 
are polynomials in E. However, the general recursion relation (23) suffers from two 
major drawbacks. In the first place, even if we choose P (E) and P^E) as polynomials 
of degree and 1 in E, respectively, (23) is incompatible with the desirable property 

that P k (E) be of degree k in E for all k, unless c = 0. Secondly, even in this case 

(23) will be in general a four-term recursion relation, implying that the polynomials 
P k (E) may not be orthogonal with respect to any (nonzero) Stieltjes measure duj(E). 
Indeed, it is well known, [2], [10], that a necessary and sufficient condition for a family 
of polynomials {Pk}^ =Q (with degP fc = k) to form an orthogonal polynomial system is 
that P k satisfies a three-term recursion relation of the form 

P k = (A k E + B k )P k _ 1 + C k P k _ 2 , k>l, (24) 

where the coefficients A k7 B k7 C k are independent of E, A k ^ 0, C 1 = 0, and C k ^ 
for k > 1. If the coefficient C k in (24) vanishes for some positive integer k, then this 
recursion relation only defines a weakly orthogonal polynomial system [11]. It is one of 
the main goals of this paper to show that both difficulties described above can always 
be overcome, provided (roughly speaking) that we expand the eigenfunction \e with 
respect to an appropriate variable. This will be achieved by using the non-uniqueness 
of -ff gauge , due to the GL(2) symmetry described in the previous section, to place H gauge 
in a suitable canonical form. 

From the form of the recursion relation (23), it follows that both difficulties de- 
scribed above disappear if 

c__ = c ++ = 0. (25) 
Indeed, if (25) holds then (23) reduces to the three-term recursion relation 



Pk-i 



-[(2k-n-\)c _ + c_]P k = 

E + c m + c (k - | - 1) + c 00 (k - ^ - l) 2 
+ (k-l)(k-2-n) [(2k - n - 3)c +0 + c + ] P k _ 2 , k > 1, (26) 

which uniquely determines all the functions P k {E) in terms of P (E) provided that, 
for all positive integer values of k, the coefficient of the left-hand side of (26) does not 
vanish. If P (E) is taken as a constant, for instance if P (E) = 1, then (26) implies 
that P k (E) is a polynomial of degree k in E for all k > 0. 

Let us see now that we can always arrange for (25) to be satisfied, by using the 
action (17) to transform P(z) into a normal form P(w) for which (25) holds. Indeed, 
(25) simply states that the polynomial P(z) vanishes at z = and z = oo, when z is 
allowed to vary over the complex projective line CP. Note that we need z to belong to 
the complex projective line at this stage so that P is guaranteed to have a root, which 
is essential for the argument that follows. Consequently, the GL(2,R) action described 
in the previous section will be replaced in what follows by a GL(2, C) action. 

We can assume, first of all, that P(z) is one of the normal forms listed in equations 
(19) and (20). We must distinguish three cases, characterized by the position of the 
roots of P in the complex projective line. Indeed, either P has two different roots 
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z 1 7^ z 2 in CP, or it has four coincident roots. In the first case, either one of the roots 
is at infinity, or both roots are finite. 

Case 1: P has two different roots z 1 7^ z 2 = 00. 
This case occurs when P is one of the first four normalizable canonical forms (19), 
or the fifth periodic canonical form (20). In this case, the translation w = z — z 1 
transforms P(z) into a polynomial P(w) vanishing at zero and infinity. In the original 
z coordinate, by (18) this amounts to replacing (21) by 

x E (z) = J2 p ^ E r J ^ lL - ( 27 ) 

In other words, we expand Xe( z ) as a P ower series around the point z = z 1 , which is a 
singular point of the linear differential equation (H gauge — E)xe = (if z i * s a simple 
root of P, z 1 is actually a regular singular point, whose indicial equation is easily seen 
to have as a root). By (7), in the "physical" coordinate x (27) becomes 

M*) = Ka*)) f2P k (E) {C{x) ~ Zl)k . (28) 

k=0 

Case 2: P has two different finite roots z 1 7^ z 2 . 
This is the case when P is one of the first four periodic normal forms (20) . The projective 
transformation w = (z — z x ) / (z — z 2 ) will again transform P(z) into a polynomial P(w) 
vanishing at w = 0, 00. Going back to the original z coordinate, by (18) we just have 
to replace (21) by 



00 , 
X E (z) = (z-z 2 rJ2y P k(E)[ 

k=0 ' ^ 



(29) 



apart from an inessential overall factor. In terms of the physical coordinate x, (29) can 
be written as 



m*) = (c(x) - z 2 r p ±p k (E) ) " ■ 



(30) 



Case 3: P has a quadruple root. 
This corresponds to the fifth normalizable canonical form, P = 1, which has a quadruple 
root at infinity. Note that P = 1 implies that the physical coordinate x can be taken 
as the canonical coordinate z. By (9), we have 

C_|__|_ = c_|_ = c 00 = c _ = 0, c — = 1. 

Performing an additional translation, if necessary, we can also take without loss of 
generality c_ = Q(0) = (notice that P is constant, and therefore does not change 
under translations). Thus equation (23) reduces in this case to 



E + c„ + c (k--) Pfc + fcCfc-l-rOc+P*.!, k>0. (31) 
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Since P = 1 is the fifth normalizable case of references [4] , [6] , c + must vanish if we want 
H to be normalizable, i.e., the algebraic eigenfunctions of H to be square-integrable. 
Therefore, in this case (31) reduces to 



~Pk+2 — 



E + c, +c (k--) ] 



Pi 



k > 0, 



which is equivalent to two two-term recursion relations for the even and odd coefficients 
Pf = P 2j and Pj = P 2j+1 , namely 



■\H-i 



E + c,+c (2j + e--) 



P e 



j > 0; e = 0, 1. 



(32) 



Note that in this case the potential is V(x) = \c\x 2 — (with c < 0), cf. [4]. 

To complete the discussion of Cases 1 and 2, we still have to deal with an important 
technical issue; namely, we must find under what conditions the coefficient of P k in (26) 
never vanishes for positive integer values of k. Let P and Q be the transforms of P 
and Q under the projective transformation z i— > w defined in the foregoing discussion of 
Cases 1 and 2; note that, by construction, P(w) vanishes at w = 0, oo. The coefficient 
of interest can be expressed as 



(2k — n — 1) c _ + c_, k > 1, 



(33) 



where 



c _ = ^P'(0), c_ = Q(0). 



From (17) it easily follows that 



1 



c _ = -P'( Zl ), 



c_ = Q(z x ) 



(34) 



for Case 1 (w = z — z x ), and 



2( Zl -z 2 y 



c = 



QM 



(35) 



for Case 2 (w = (z — z 1 )/(z — z 2 )). We shall now distinguish three subcases: 

Case i. z 1 is a simple real root of P 
This case occurs when P is one of the canonical forms 2, 4, 6, 7, or 10. Note that in 
this case the mapping z i— > w is real, and so are the coefficients c _, c_. From (12) and 
(34)-(35) it is immediate to deduce the asymptotic formulas 



x ~ \z 



1 

1 2 



H(z) 



+ 21 



where we have dropped unessential constant multiplicative factors from the right-hand 
side, and have taken for convenience z x as the lower limit of the integral giving x 
in terms of z. We saw in the previous section that when H is fully algebraic it has 
n + l linearly independent algebraic eigenfunctions of the form (7), where x <^V n . It 
follows that the polynomial factor x(z) cannot vanish at the origin for all the algebraic 
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eigenfunctions of H. Hence there is at least one algebraic eigenfunction of H whose 
asymptotic behavior at x = is given by 



*(£r-»). 

If all the algebraic eigenfunctions of H are regular at x = 0, then we must have 



ip(x) ~ \x\ 



-o- 



n > 0. 



(36) 



Since (33) can be written as 



2Cn 



n + 



k- 1 - 

2 



it follows from (36) that the coefficient (33) cannot vanish in this case. 

Case ii. z 1 is a simple complex root of P 
In this case P is either the first or the eighth canonical form. Since z 1 is not real, the 
mapping z i— > w is not real either, and the above asymptotic argument is not valid (the 
eigenfunctions of H need not be regular outside the real axis). For the first canonical 
form (19), we can take w = z — i and therefore 



i v, 



- c + + ic 



from (34). Hence the coefficient (33) does not vanish in this case provided that the 
following conditions are satisfied: 



or 



i( n+ l-^l,2,.... 



(37) 



It is easily checked that the choice w = z + i leads exactly to the same conditions. For 
the eighth canonical form, we can take w = (z — i)/(z + i), and therefore, from (35), 



1 2 
-UK 

2 



= ^ + i( c ,-c ). 

2 2 V + 



Hence in this case the conditions for the coefficient (33) not to vanish are 
c_^c + or I( n+ 1--£^1,2,.... 



(38) 



It is straightforward to check that the choice w = (z + i)/(z — i) yields the same 
conditions, while the other natural choice w = (y/l — k 2 z =f i)/ (Vl — k 2 z ± i) only has 
the effect of replacing the first condition (38) by c + ^ (1 — k 2 )c_. 

Case Hi. z 1 is a multiple root of P 
This case takes place when P is either the third or the ninth canonical form, and in 
both cases (33) reduces to c_. For the third canonical form (19), if c_ 7^ then we 
take w = z, and therefore c_ = c_ 7^ 0. If c_ = 0, then c + 7^ if all the algebraic 
eigenfunctions of H are square-integrable (see [4]). Hence, taking w = 1/z, we get 
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P(w) = vw 2 and Q = — (c + + c w), so that c _ = and c_ = — c + 7^ 0. Hence the 
coefficient (33) cannot vanish in this case. Finally, if P is the ninth canonical form (20) 
then w = (z — i)/(z + i) and 



C- = C i + kc + -c_). 



2 2 

Hence (33) will not vanish if 

c + 7^ c_ or c 7^ 0. 
Note that when (39) does not hold V reduces to a constant potential: 



(39) 



v = T v -T2 n(n + 2) - c - 

In summary, the previous analysis shows that the critical coefficient (33) cannot 
vanish for any positive integer k provided that V is fully algebraic, that all its algebraic 
eigenf unctions are regular (or square-integrable, for the third normalizable canonical 
form (19)), and that conditions (37), (38), and (39) are satisfied when P is one of the 
normal forms 1, 8 or 9, respectively. If (33) doesn't vanish, defining new polynomials 

(-i) fc h 



, J (26 °- )fe r(^ + *-i + |) 
' (-i) fc s 



if c _ = 



the recursion relation (26) can be written in the more standard form 



(40) 



fc+i 



E + c^+c (fc-^)+c 00 (k-^)' 



-k(k-n-l)[c +0 (2k-n-l) + c + ][c _(2k-n-l) + c_]P k _ 1 , k>0. (41) 
We have thus proved the main theorem in this section: 

Theorem 1. Let V be a fully algebraic one-dimensional quasi-exactly solvable poten- 
tial whose algebraic eigenfunctions are all regular (or normalizable, if V corresponds 
to the third or hfth canonical forms in (19)). Assume, furthermore, that conditions 
(37), (38) or (39) are satisfied, if V is obtained from the first, eighth or ninth canonical 
forms (19)-(20), respectively. Then V defines a family of weakly orthogonal polynomi- 
als {-PfcjfcLo sa ^ s fy m S a three-term recursion relation (41) (or (32), if V corresponds 
to the hfth canonical form). The polynomials P k are defined by (40) and (28), if V 
is associated to one of the canonical forms 1-4 or 10, or by (40) and (30), if V corre- 
sponds to one of the normal forms 6-9. Finally, the potential V associated to the hfth 
canonical form defines two families of weakly orthogonal polynomials P® = P 2 j and 



P 



P 2j+1 through (28). 
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We shall say that a quasi-exactly solvable potential V is exactly solvable if it is 
independent of the "spin" parameter n. This implies that V has n algebraic eigenval- 
ues and eigenfunctions for arbitrary n G N, so that we can algebraically compute an 
infinite number of eigenvalues of V (leaving aside the boundary conditions). All ex- 
actly solvable normalizable one-dimensional potentials have been classified; see [4] for 
a complete list. The quintessential example of exactly solvable one-dimensional poten- 
tial is the harmonic oscillator potential, which corresponds to the fifth canonical form 
(19). We have seen in the previous section that in this case there are two families of 
orthogonal polynomials (the odd and even coefficients in (28)), each of which satisfies 
a two-term recursion relation (32). We shall now show that, as conjectured in [1], the 
latter property actually characterizes exactly solvable normalizable potentials: 

Theorem 2. The weakly orthogonal polynomial system associated to an exactly solv- 
able normalizable potential satishes a two-term recursion relation. 

Proof: The proof is a simple case- by-case analysis using the classification of exactly 
solvable normalizable potentials given in [4]. Indeed, for the first normalizable canonical 
form P(z) = u(z 2 + 1) we have w = z^fi, and therefore P(w) = P(w±i) = vw(w±2i), 
so that c +0 = 0. Since Q(w) = Q(w±i), we also have c + = Q"(0)/2 = Q"(±i)/2 = c + . 
But the exactly solvable potentials associated to this normal form are characterized 
by the vanishing of c + , [4], so that c + = c +0 = 0, and (41) is a two-term recursion 
relation. Similarly, for the second normalizable canonical form, P(z) = u(z 2 — 1) and, 
for instance, w = z =F 1- Proceeding as before we obtain that c +0 = and c + = c + . 
Since exactly solvable potentials are again those satisfying the condition c + = 0, (41) 
reduces to a two-term recursion relation. 

The third normalizable canonical form has P(z) = vz 1 , and therefore c +0 = c _ = 
0. The exactly solvable potentials are characterized by the vanishing of the coefficients 
c + or c_, but not both simultaneously. In the former case we can take w = z, while 
in the latter w is proportional to 1/z (see the foregoing discussion on the vanishing 
of the critical coefficient (33)). In either case, the coefficient of P k _x in (41) vanishes 
identically. 

The fourth normalizable canonical form is given by P(z) = z, so that w = z and 
c_|_ = c +0 = 0, and its exactly solvable potentials are defined by the vanishing of the 
coefficient c + = c + = 0, so that (41) is two-term. Finally, for the fifth normalizable 
canonical form P(z) = 1 all normalizable potentials are automatically exactly solvable 
(they are translates of the harmonic oscillator), and we have already seen that its 
associated orthogonal polynomials satisfy the two-term recursion relations (32). Q.E.D. 

4. The Orthogonal Polynomials 

We shall study in this section the properties of the family of weakly orthogonal 
polynomials associated to a quasi-exactly solvable one-dimensional Hamiltonian in the 
manner described in the previous section. Since, as we shall see, these properties can 
be established directly from the recursion relation (41) or (32), these polynomials have 
basically the same properties as those studied by Bender and Dunne in [1] . 

We have seen in the previous section that the polynomials P(E) defined by a quasi- 
exactly solvable one-dimensional Hamiltonian satisfy a three-term recursion relation of 
the form 

n + i = ( £ -Wn-%^-i. k >^ ( 4 2) 
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with a = and 

a n+1 = 0. (43) 

For the fifth canonical form, the polynomials P k and P k also satisfy a recursion relation 
of the form (42), with a k = for all k > 0. Note that the coefficients a k , b k in (42) are 
guaranteed to be real only for the canonical forms 2-7 and 10 (for which P has a real 
root). As remarked in the previous section, the vanishing of a k for a positive integer 
value of k means that the polynomials P k are only weakly orthogonal. In particular, 
many classical results, based on the fact that a k > (or sometimes a k > 0) for k > 1 
cannot be applied in our case. 

By Favard's theorem, [2], there is a moment functional, that is a linear functional C 
acting in the space C[E] of (complex) univariate polynomials, such that the polynomials 
P k are orthogonal under C: 

£(P k P l ) = lk 5 kl , kjEN. (44) 

The functional C is unique if we impose the normalization condition C(P ) = = 1. 
It is also known (Boas's theorem, [2]) that there is a (not necessarily unique) function 
of bounded variation u such that 

/oo 
p(E)du(E) (45) 
-oo 

for an arbitrary polynomial p. The coefficient 7 fc = C{P 2 k ), which therefore plays the 
role of the square of the norm of P k , can be computed by multiplying (42) by P k _ 1 
and taking C of both sides, obtaining 



° = 7fc-a fc 7fc-i, k>l. 
Taking into account that 7 = C(l) = 1 we get 

k 

lk = U a P k > 1. (46) 

In particular, from this formula follows one of the key properties of the weakly orthog- 
onal polynomial system associated to a one-dimensional quasi-exactly solvable Hamil- 
tonian. Namely, from (43) we have 

7fc=0, k>n+l, 

so that all the polynomials P k with k > n + 1 have zero norm. From this formula it also 
follows that the "squared norms" 7 fc will be positive for k < n if and only if a k > for 
1 < k < n. It can be shown by a straightforward computation that this is always the 
case when P is one of canonical forms 2-4 in (19), assuming that all the eigenfunctions 
of H are square-integrable and that H is not exactly solvable. Note also that when H 
is normalizable (canonical forms 1-5 in (19)) and exactly solvable then a k = for all 
k > 0. Hence the square norms of all the polynomials P k vanish, from which it easily 
follows from (42) that £ = S(E - b ). 
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Other important properties of the polynomials P k concern their zeros. Classically, 
[2], it can be shown that if a k > for all k e N then the zeros of the polynomials 
P k satisfying a three-term recursion relation (42) are real and simple. In our case the 
condition a k > for all k G N can never hold on account of (43). However, if H is fully 
algebraic it can still be proved that all the zeros of P n+1 are real and simple. Indeed, by 
hypothesis H is self-adjoint on the space 9T of functions of the form (7), with x £ P n - 
Hence H has n + 1 linearly independent algebraic eigenf unctions lying in 9t, whose 
corresponding eigenvalues are real (by self-adjointness) and distinct (H being a one- 
dimensional Sturm-Liouville operator) . Let us denote by E < E x < . . . < E n these 
n+1 real eigenvalues of H on 91, and by i(^i(x) = ifj Ei (x) the eigenfunction corresponding 
to the eigenvalue E x . Then (7) and either (28) or (30) imply that P k (E t ) = 0, or 
equivalently P k {E{) = 0, for k > n + 1 and < / < n. In particular, since P n+1 is of 
degree n + 1 and all the eigenvalues E l are different, it follows that 

n 

P n+1 (E) = H(E-E l ), (47) 

1=0 

where we have used (42) and the fact that P Q = 1. In other words, P n+1 has n + 1 
simple real zeros at the n + 1 algebraic eigenvalues of H . Furthermore, from the fact 
that P k vanishes at E x for k > n + 1 we conclude that there exist monic polynomials 
Q k of degree k such that 

P k+n+1 =Q k P n+1 , k>0. (48) 

This is the so called factorization property of the polynomial system {-P/jjfcg^' °f- [!]• 
Note that the vanishing of P k (E t ) for all k > n + 1 is consistent with the recursion 
relation on account of (43). In fact, when a k is positive for k > 1 and b k is real for k > 0, 
(47) follows directly from the recursion relation by Lemma 3, without using the fact 
that the polynomials P k are associated to a fully algebraic quasi-exactly solvable one- 
dimensional Hamiltonian. The vanishing of P fc (-E ; ) for k > n + 1 is then an immediate 
consequence of P n+1 (E' Z ) = 0, the recursion relation (42) and (43). 

From the previous equation and (42) it follows that the polynomials Q k also satisfy 
a three-term recursion, namely 

Q k+1 = (E- b k+n+1 )Q k - a k+n+1 Q k _ l7 k>0, 

and are therefore orthogonal with respect to an appropriate moment functional Cq (in 
general different from C). 

It was heuristically argued in [5] that 

n 

C = Y t u j S{E-E j ) (49) 

3=0 

on C[E], where the coefficients uj- are defined by 

n 

J2P k (Ei)^ = 5 k0 , fc = 0,l,...,n. (50) 

1=0 
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Equivalently, the discrete Stieltjes measure du(E) defined by the function 

n 

u>(E) = J2" j 8(E-E j ), (51) 

j=0 

where 6{t) is Heaviside's step function, satisfies (45). Note that the linear system (50) 
uniquely defines the n + 1 constants ujj, since by (27) or (29) its coefficient matrix is the 
matrix of the change of basis {c k (z — z x ) k /k\} k _ Q or {c k (z — z 1 ) h (z — z 2 ) n ~ k / k _ Q 

to {xei^-q * n ^-'®^ > ni c k being the coefficient of P k in (40). It is not difficult to show 
rigorously that (44) is satisfied. Indeed, by the uniqueness of C this is equivalent to 
showing that if C = Ej=o u j — E-) then 

£ o (P fe P,) = 0, Ml, (52) 

and that 

£ (P )=£ (1) = 1, 

since C{P\) and £ (P|) must coincide if (52) holds due to the recursion relation (42). 
From the definition of uj- we deduce that the last equation, together with (52) for 
k = and / = 1, . . . , n, are satisfied. Suppose now that (52) holds for k = 0, 1, . . . , K 
[K < n — 1) and k < I < n. Multiplying (42) by P l and taking C of both sides we 
obtain 

£ (P K+1 P Z ) = C ((E - b^PxPi) - a K £ (P i ,_ 1 P z ) = Co(EP K Pi) 
if K + 1 < / < n, by the induction hypothesis. But, using again (42), 

£ (EP K Pi) = £o(Pk ■ EPi) = C (P K P l+1 ) + ^(PtfP,) + a^pyVj = 0, 

by the induction hypothesis (since I > K + 1 implies / — 1 > K). Hence (52) is true 
for < k, I < n. Finally, (52) is trivially true when k or / are greater than n by the 
factorization property (48) and (47) . 

We shall next show that all the coefficients Uj are positive if b k is real for all 
< k < n and a k > for 1 < k < n. (Several instances of this result were checked 
numerically in [5] for the orthogonal polynomials associated to the Hamiltonian (1).) 
The proof is based on the following simple lemma: 

Lemma 3. If a k > for k = 1, 2, ... ,n and b k is real for k = 0,1, ... ,n then C is 
positive-definite on V 2n - In other words, if p G V 2n is a real polynomial of degree at 
most 2n,p^0 and p(E) > for all E e R then C(p) > 0. 

Proof: A polynomial p G V 2n which is non-negative for all real values of E must be 
of the form q 2 +r 2 , where q,r E V n are real polynomials. Write q = ^2 k=0 %Pk', then all 
the coefficients q k are real, since P k is a real polynomial for < k < n by the hypotheses. 
Using the orthogonality of the polynomials P k we obtain C(q 2 ) = Y^k=o Qk'fk- Similarly, 

if r = Y! k =oPkPk then C ( r ) = ELo r fc7 fc , and C(p) = ELo(^ + r lhk- Since 7 fc > o 
for k = 0, 1, . . . , n by (46) and the hypothesis on the coefficients a k , it follows that 
C(p) > 0, and C(p) = if and only if q k = r k = for k = 0,1,..., n, that is if 
p = 0. Q.E.D. 
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Proposition 4. If a k > for k = 1, 2, . . . , n and b k is real for k = 0, 1, . . . , n then 
uj k > for all k = 0, 1, . . . n. 

Proof: Apply the previous lemma to the polynomials Y\_ Q< j^ k<n {E ~~ -^j) 2 e ^2n 
for k = 0, l,...n. Q.E.D. 

Note that the hypotheses of the previous proposition are satisfied when P is one of 
canonical forms 2, 3 or 4, provided that all the eigenfunctions of H are square-integrable 
and that H is not exactly solvable. In particular, it is satisfied by the Hamiltonian (1). 

The (Hamburger) moment problem for the moment functional (49) associated to 
the weakly orthogonal polynomials defined by a quasi-exactly solvable one-dimensional 
Hamiltonian consists in determining whether there is a distribution function (i.e., a non- 
decreasing function of bounded variation) uj such that C can be represented by (45) for 
an arbitrary polynomial p. We have already shown that this problem has a solution 
(51), since (51) is clearly non-decreasing and of bounded variation. We shall next show 
that this solution is unique (up to an additive constant), so that the moment problem 
associated to the weakly orthogonal polynomial system {-P^jyteN ^ s always determined 
[12]. Essentially, this is due to the fact that the spectrum 

a(u) = {EeR: uj(E + S) - uj(E - S) > 0, V5 > 0} 

of the distribution function (51) is the finite set {Ej}™ =G [13]. According to a well known 
result in the classical theory of orthogonal polynomials, [2], a distribution function uj 
defines a positive-definite functional on C[E] through integration with respect to the 
Stieltjes measure duj(E) if and only if the spectrum of uj is infinite. Since C is not 
positive-definite (£(P^ +1 ) = 7 n+1 = 0), any solution uj of (45) must have a finite 
spectrum, and will thus be of the form 



uj(E) = J2" k e(E-E k )+C 



k=0 



for some constant C, up to an immaterial redefinition of uj in a(uj). If / is a compact 
interval containing cr(cD) U a(u), then 

C(p) = Jp{E) du{E) = Jp{E) dw(E), Vp G C[E]. 

Since / is compact, a well known theorem (cf. [2]) shows that uj and u differ by a 
constant at all points in which both uj and u are continuous. But this easily implies 
that E k = E k and Lo k = Hi k for k = 0, 1, . . . , n = n, whence uj = uj + C, as stated. Note 
that the same argument shows that the moment problem in any interval containing 
[E Q ,E n ]; in particular, the (Stieltjes) moment problem in [E Q ,oo) is also determined. 
In this respect, the weakly orthogonal polynomials associated to a quasi-exactly solvable 
one-dimensional Hamiltonian behave in exactly the same way as the classical orthogonal 
polynomials, whose moment problem is also determined, [2]. 

The moments of the moment functional C are by definition the numbers 



/oo n 
E k duj(E) = J2 u i E t k eN. (53) 
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If the hypotheses of Proposition 4 hold, all the moments are real. From (53) we see 
that the module of the k-th moment \i k does not grow factorially as k tends to infinity, 
as argued in [1], but instead it diverges like the k-th power of a constant [14]. 
We shall next show that if the coefficient a k satisfies the condition 

a k ^ 0, 1 < k < n, (54) 

which guarantees that the polynomials P k have non-zero norm for k < n, then the 
moments \x k with k > n + 1 satisfy a constant coefficient difference equation of order 
n + 1. To this end, recall first of all that the bilinear form (p,q) = C(pq) defined 
by C in C[E], when restricted to the subspace C <g> 7^, is represented in the basis 
{£' fc } 0<fc< ; by the symmetric matrix (fi i+ j) 0<i J<z , whose determinant we shall denote 

by A;. On the other hand, the matrix of the bilinear form (• , •) in the basis {Pk}o<k<i 
is clearly diag(l, 7 l5 . . . , 7 ; ); therefore, by (46) and the hypothesis on the coefficients 
a fc , we conclude that A n ^ and 

A fc = 0, k>n + l. (55) 

In particular, since A n ^ but A n+1 = 0, the last column of A n+1 must be a linear 
combination of the remaining columns, so that 

n+1 

i"fc = I>^-*' n + l<k<2(n + l), (56) 

i=i 

for some (in general complex) constants c 1 , . . . , c n+1 . An easy induction argument 
using (55) then shows that the above relation is actually valid with the same constant 
coefficients c i for all k > n + 1, as claimed. In fact, it is not hard to see that c i in (56) 
is minus the coefficient of E n+1 ~ % in P n+1 . Indeed, write P n+1 = E n+1 — p n , with 

n+1 

Pn = E C ^ n+1 " 1 ' 
i=l 

and let Q k = E k — q k _ x , so that q_ x = and degq k _± < k — 1 for k > 1. From (48) it 
follows that 

P k = E k - E k - n ~ l p n - q k _ n _ 2 P n+1 , k>n+l, 
which by (44) implies that 

n+1 

^ k = C(E k ) = CiE^-'pJ = J2c t ^ n + l<k<2(n + l). 

i=i 

Comparing with (56) and taking into account the linear independence of the columns 
of A n we immediately obtain that c i = c i for % = 1, 2, . . . , n + 1, as stated. 

Note that the fact that the moments satisfy a constant coefficient recursion relation 
(56) (with k > n + 1) actually characterizes weakly orthogonal polynomial systems. 
Indeed, (56) simply expresses the fact that the (n + 2)-th column of A ; for / > n + 1 
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is a linear combination of the first n+1 columns. Hence the recursion relation implies 
(55), and since A n+1 = n?=i 7? this m eans that "f k = for some k < n + 1, so that 
a fc = for some /c < n + 1 by (46) . 

Consider, for example, the Hamiltonian (1) studied in [1], which corresponds to 
the fourth canonical form with 

n = J — 1, c + = —16, c = = 0, c_ = 2s + -(n — 1). (57) 

The coefficients of the corresponding recursion relation (42) are easily found to be 

6 fc = 0, a k = 16k(J -k)(k + 2s- 1), fc > 0. (58) 

Since we can take s > 1/2 without loss of generality, we see that a fc > for 1 < k < n, 
so that (54) is satisfied. Furthermore, since b k vanishes for all k > the polynomials 
P k have parity (— l) fc , and therefore all the odd moments vanish (the corresponding 
moment functional is said to be symmetric). For J = 3 (that is, n = 2), according to 
the foregoing observations we know that the moments satisfy a third-order recursion 
relation of the form (56), whose coefficients are minus the coefficients of E 2 , E and 1 
in P 3 . From (42) (with P Q = 1) we obtain 

P X =E, P 2 =E 2 -64s, P 3 (E) = E 3 -32(4s + l)E, (59) 

so that c 1 = c 3 = — as expected, since the moment functional is symmetric — , and 
c 2 = 32(4s + 1). Therefore the even moments satisfy the first-order recursion relation 

/ u 2j =32(4s + 1) / u 2j _ 2 , j>2, (60) 

and since (i 2 = r y 1 = a x = 64s, from (60) we obtain 

fi 2j = 32 J ~ 1 (4s + ly- 1 ■ 64s, j > 1. (61) 

Thus, in this case fx 2 j has a pure power growth. The same result can be obtained using 
(53). Indeed, from (59) we have 

E = -A = -^32(4^ + 1), E 1 =0, E 2 = A, 

and therefore 



4s + 1 4s + 1 1 



s 2s + 1 

= Is 

from (50) and (59). Thus 

= 47TI [{ ~ x)k + xk] 

which yields /U 2j+1 = for j > and (61). 
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5. Conclusions 



We have shown in this paper how every quasi-exactly solvable one-dimensional 
Hamiltonian satisfying conditions (37)-(39) defines a weakly orthogonal polynomial 
system {P k }kL through the three-term recursion relation (41) (with initial condition 
Pq = 1). It is important, in this context, to emphasize the weak orthogonality of the 
polynomials P k , i.e., the fact that the norm of P k may vanish — and in fact does vanish 
for k > n + 1, n being the "spin" parameter present in the Hamiltonian. As explained 
in Section 4, this is an inevitable consequence of the vanishing of the coefficient of 
P k _ 1 in the recursion relation (41) for k = n + 1, which is made possible by the fact 
that the parameter n is a non-negative integer. The latter fact, however, is an intrinsic 
property of one-dimensional quasi-exactly solvable (as opposed to merely Lie-algebraic) 
Hamiltonians; indeed, it is a key factor in the explanation of the partial integrability of a 
quasi-exactly solvable Hamiltonian outlined in Section 2. To better illustrate this point, 
consider the Hamiltonian (1), which is Lie-algebraic for all real values of the parameter 
J. Indeed, H can be written in the form (5)-(6), with ((x) = x 2 /4, fj,(z) = e~ 4z z s_1 / 4 , 

c ++ = c +0 = c 00 = c H = c__ = 0, c _ = 1/2, and the remaining coefficients given 

by (57), where now n is to be regarded as an arbitrary real parameter. When n is not 
a non-negative integer, the generators (4) don't leave invariant any finite-dimensional 
polynomial module V n , so that H is in general non-integrable — there is no special 
reason for H to have algebraically computable eigenfunctions of the form (7), with x 
a polynomial. However, even when n is not a non-negative integer, the Lie-algebraic 
nature of H and conditions (25) imply that the polynomials P k defined by (7), (21) and 
(40) still satisfy a three-term recursion relation (42), with the coefficients given by (58). 
In other words, what makes H quasi-exactly solvable is not merely the fact that its 
associated polynomials satisfy a three-term recursion relation (42) (which implies their 
orthogonality with respect to some Stieltjes measure) , but the fact that the coefficient 
a k in this recursion relation vanishes for some positive integer value of /c, so that the 
associated polynomials P k can only be weakly orthogonal. 

As we saw in Section 4, the Stieltjes measure with respect to which the polynomials 
P k associated to a quasi-exactly solvable Hamiltonian H are orthogonal is supported in 
the set of algebraic eigenvalues of H, which is a finite set. For this reason, the polynomi- 
als P k are discrete polynomials. Although the classical (Hermite, Legendre, Laguerre, 
Tchebycheff, etc.) polynomials of Mathematical Physics are orthogonal with respect 
to a continuous measure, discrete (Charlier, Hahn, Krawtchouk, Meixner, Tchebycheff, 
etc.) polynomials have also been studied in the mathematical literature of orthogonal 
polynomials, cf. [2]. Note that a discrete polynomial system is truly — as opposed to 
weakly — orthogonal if and only if the supporting set of its Stieltjes measure is infinite. 
Some of the discrete polynomials cited above, like the Hahn, Krawtchouk or discrete 
Tchebycheff polynomials, are in fact weakly orthogonal. In general, weakly orthogonal 
polynomials arise naturally, for instance, in the theory of approximate polynomial curve 
fitting, [15]. More recently, [16], the study of second-order finite difference eigenvalue 
equations with infinitely many polynomial solutions has led to an interesting connection 
between a non-standard finite-dimensional representation of sl(2) and certain families of 
weakly orthogonal discrete polynomials (Hahn polynomials and analytically continued 
Hahn polynomials). 

Let us stress, in closing, that the present paper deals only with one-dimensional 
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quasi-exactly solvable Hamiltonians. It is an interesting open problem to generalize 
these results to quasi-exactly solvable multi-dimensional systems, a possibility already 
considered in [5] , where a heuristic (but inconclusive, in our opinion) argument was ad- 
vanced suggesting that all quasi-exactly solvable systems give rise to weakly orthogonal 
polynomials. In the two-dimensional case, at least, the classification of quasi-exactly 
solvable Lie algebras of first-order differential operators in two variables presented in 
[7] and [17] could be used as a starting point for an analysis along the present lines. 
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